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ABSTRACT 

We present Spitzer/IRS spectral mapping observations of the luminous infrared galaxy (LIRG) 
Arp 299 (IC 694 + NGC 3690) covering the central ~ 45" ~ 9 kpc. The integrated mid-IR spectrum 
of Arp 299 is similar to that of local starbursts despite its strongly interacting nature and high infrared 
luminosity, Ltr ~ 6 x 10 11 Lq. This is explained because the star formation (probed by e.g. high 
[Nelll] 15.56 fim/ [Nell] 12.81 /mi line ratios) is spread across at least 6-8 kpc. Moreover, a large fraction 
of this star formation is taking place in young regions of moderate mid-IR optical depths such as the 
C+C complex in the overlap region between the two galaxies and in H n regions in the disks of the 
galaxies. It is only source A, the nuclear region of IC 694, that shows the typical mid-IR characteristics 
of ultraluminous infrared galaxies (ULIRGs, Lir > 1O 12 L ), that is, very compact (less than lkpc) 
and dust-enshrouded star formation resulting in a deep silicate feature and moderate equivalent widths 
of the PAHs. The nuclear region of NGC 3690, known as source Bl, hosts a low-luminosity AGN and 
is surrounded by regions of star formation. Although the high excitation [NeV] 14.32 /jm line typical 
of AGN is not detected in Bl, its upper limit is consistent with the value expected from the X-ray 
luminosity. The AGN emission is detected in the form of a strong hot dust component that accounts 
for 80 — 90% of the 6 jiva. luminosity of Bl. The similarity between the Arp 299 integrated mid-IR 
spectrum and those of high-z ULIRGs suggests that Arp 299 may represent a local example, albeit 
with lower IR luminosity and possibly higher metallicity, of the star-formation processes occurring at 
high-z. 

Subject headings: galaxies: evolution — galaxies: nuclei — galaxies: Seyfert — galaxies: structure — 
infrared: galaxies — galaxies: individual: Arp299 
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1. INTRODUCTION 



ill 



Luminous Infrared (IR) Galaxies (LIRGs, Lm = 10 1 
10 12 Lq, see Sanders & Mirabel 1996 for a review) are an 
important cosmological class as they are the main con- 
tributors to the co-moving star formation rate density of 
the universe at z ~ 1 (Elbaz et al. 2002; Le Floc'h et al. 
2005; Perez-Gonzalez et al. 2005; Caputi et al. 2007). 
At higher redshifts (z ~ 2) LIRGs and Ultra Luminous 
Infrared Galaxies (ULIRGs, L m = 10 12 - 10 13 L Q ) have 
comparable contributions. Moreover, a large fraction of 
massive galaxies at high redshift are LIRGs and ULIRGs, 
indicating that a phase of intense star formation and/or 
AGN activity dominated the energy output of these mas- 
sive galaxies and was more common in the past than now 
(Caputi et al. 2006). 

Recent results with the Infrared Spectrograph (IRS) 
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instrument (Houck et al. 2004) on Spitzer are now show- 
ing that the mid-IR spectra of z ~ 1.7 — 2 star-forming 
galaxies with high IR luminosities (in the ULIRG class, 
or even more luminous) appear to be more consistent 
with those of local starbursts and LIRGs, rather than 
local ULIRGs (Rigby et al. 2008; Farrah et al. 2008). 
In particular these high-z IR galaxies have similar PAH 
features to those of local ULIRGs, but substantially shal- 
lower 9.7 /xm absorption features. This may be due to 
lower optical depths, perhaps because star formation in 
ULIRGs was more spatially extended at high-z than now, 
and/or metallicity effects. Therefore, local LIRGs play a 
fundamental role for understanding their more distant 
and luminous counterparts. To characterize fully the 
mid-IR properties of local LIRGs, we have started a pro- 
gram to obtain IRS spectral mapping of a sample of lo- 
cal LIRGs from Alonso-Herrero et al. (2006). In this 
paper we present the first results for Arp 299 (IC 694 + 
NGC 3690), a LIRG with L m ~ 6 x 10 n L Q , using the 
IRAS fluxes from Sanders et al. (2003) and a distance of 
d = 42Mpc (Ho = 75kms _1 Mpc _1 ). 

Arp 299 has long been recognized as an extremely in- 
teresting system: an interacting galaxy (Arp 1966); a 
non-Seyfert 9 galaxy with bright emission lines (Weed- 
man 1972); a luminous mid-IR source (Rieke & Low 
1972); and a highly powerful starburst (Gehrz, Sramek, 
& Weedman 1983). The bright sources in this system 
probe a full range of distinct and well-characterized phys- 

9 but see Garci'a-Marm et al. 2006 for the optical identification 
of the AGN in NGC 3690 
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TABLE 1 

Log of the Spitzer/IRS Spectral Mapping Observations 



Target 


RA 


Dec 


IRS 


Map 


Plate 


PA 


Ramp 


Number 


Program 








Module 


Size 


Scale 




Duration 


Cycles 


ID 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Arp 299 


Ilh28m32.8s 


+58d33m42s 


SL1 + SL2 


42 x 19 


1.85 


56.82 


14s 


2 


21 








LL1 + LL2 


37 x 8 


5.08 


-26.9 


14s 


2 


21 


IC 694 


Ilh28m33.7s 


+58d33m46s 


SH 


16 x 10 


2.26 


157.38 


30s 


2 


30577 








LH 


5x6 


4.46 


72.58 


60s 


4 


30577 


NGC 3690 


Ilh28m30.7s 


+58d33m43s 


SH 


16 x 10 


2.26 


157.31 


30s 


2 


30577 








LH 


5x6 


4.46 


72.51 


60s 


4 


30577 



Note. — (1) Target; (2) and (3) RA Dec (J2000) of the center of the final data cube (before rotation) as given in the headers; (4) IRS 
Module; (5) Original size of the final map in pixel before rotation; (6) Plate scale of the map in arcsec/pixel; (7) Average slit orientation 
over the DCEs in degrees; (8) Ramp duration; (9) Number of cycles; (10) Spitzer archive Program ID. 



ical conditions that can be used as a benchmark for dif- 
ferent mid-IR indicators for LIRGs. In brief, the prop- 
erties of the nuclear region of IC 694, referred to as A 
(see Gehrz et al. 1983), are dominated by a relatively 
extended and evolved period of star formation, whereas 
the C+C complex in the overlap region between the two 
galaxies, contains powerful and young starbursts. The 
nuclear region of NGC 3690, referred to as Bl, hosts an 
obscured AGN and is also surrounded by regions of star 
formation. §5, 6, and 7 discuss the properties of these 
sources. Prior to that discussion, §2 presents the ob- 
servations and data reduction, whereas §3 presents the 
analysis of the data. In §4 we describe the overall mor- 
phological mid-IR properties of the system. Finally, §8 
puts the mid-IR properties of Arp 299 in the context 
of local and high-z IR-bright galaxies, and §9 gives the 
conclusions. 

2. OBSERVATIONS 
2.1. IRS Spectral Mapping Observations 

The observations of the interacting system Arp 299 
were obtained with IRS using all four modules: Short- 
High (SH; 9.9 - 19.6 ^m), Long-High (LH; 18.7 - 
37.2 urn), Short-Low (SL1; 7.4 - 14.5 /im and SL2; 5.2 - 
7.7/xm) and Long-Low (LL1; 19.5 — 38 /im and LL2; 
14.0 — 21.3 /um). The low resolution modules contain two 
subslits, one for the first order (SL1 and LL1) and the 
other for the second order (SL2 and LL2). The high res- 
olution modules are cross dispersed so that ten orders 
fall on the array (see Houck et al. 2004 for full details). 
The IRS high (R ~ 600) spectral resolution data (SH 
and LH modules) were obtained in Cycle 3 as part of a 
larger GTO program (P.I.: G. H. Rieke) aimed at ob- 
taining Spitzer/IRS spectral mapping of a representative 
number of local LIRGs from the volume-limited sample 
described by Alonso-Herrero et al. (2006). The spectral 
mapping data in low (R ~ 60 — 126) spectral resolution 
(SL and LL modules) were taken from the Spitzer archive 
and were part of a different GTO program (P.I.: J. R. 
Houck) observed in Cycle 1. 

The observing strategy of both programs used the IRS 
spectral mapping capability. In this observing mode the 
telescope is moved perpendicular to the long axis of the 
slit with a step of one-half the slit width until the appro- 
priate region is covered. This strategy produced a map 
redundancy of 2 for the SL, LL, and LH modules, and of 4 
for the SH module. The two galaxies of Arp 299 were ob- 
served as separate targets with the SH and LH modules, 
whereas in low resolution they were observed with one 



pointing. Table 1 gives details of the IRS spectral map- 
ping observations, including the ramp integration times, 
number of cycles, and final sizes of the data cubes. 

Since the SL and LL slits arc longer than the extent 
of the galaxy, no separate background observations were 
needed. For the LH module we obtained dedicated back- 
ground observations with a single slit (staring mode) ob- 
servation at a region about 2 arcminutes away from the 
galaxy. No backgrounds were observed for the SH mod- 
ule. We found however, that even for the LH module 
where the background is higher than for the SH module, 
the nuclei are so bright that the background subtraction 
had little effect. 

The data were processed using the Spitzer IRS pipeline 
version SI 5. 3. The typical uncertainties of the wave- 
length calibration arc 0.009 /im for SL1, 0.006 /mi for 
SL2, 0.036 /im for LL1, 0.034 /im for LL2, 0.003 /mi for 
SH, and 0.01 /im for LH. The uncertainty of the absolute 
photometric calibration is less than 10%. The observa- 
tions are diffraction limited above 6 /im, so the spatial 
resolution depends on the wavelength of the observation 
and it is generally 1.5 — 2.5 pixels/spatial resolution el- 
ement (see IRS Handbook, version 3.1). Our own es- 
timates of the angular resolution (FWHM) of the data 
cubes are typically 4" for SL2, 5" for SL1, and 5.5" for 
SH (see Pereira-Santaella et al. 2009, in preparation). 
For the assumed distance to Arp 299 these correspond 
to angular resolutios of 0.8 kpc, lkpc, and 1.1 kpc for 
the SL2, SL1, and SH spectral maps, respectively. 

The data cubes were assembled using cubism (the 
CUbe Builder for IRS Spectra Maps, Smith et al. 2007b) 
from the individual Basic Calibrated Data (BCD) spec- 
tral images. Full error cubes are also built alongside the 
data cubes by standard error propagation, using, for the 
input uncertainty, the BCD-level uncertainty estimates 
produced by the IRS pipeline from deviations of the fit- 
ted ramp slope fits for each pixel. These uncertainties 
are used to provide error estimates for extracted spectra, 
and constructed line and continuum maps (see Smith et 
al. 2007b for full details). These errors together with the 
additional systematic uncertainties arising from calibra- 
tion effects and pointing errors bring the error budget to 
typically 15 - 20%. 

2.2. IRS Staring Spectroscopy 

Additional IRS staring mode spectroscopy with the 
SH and LH modules of the two nuclei of Arp 299 was 
obtained from the Spitzer archive. These observations 
were observed as part of Program 21, and had no as- 




arcsec arcsec 

Fig. 1. — Observed (not corrected for extinction) SL module spectral maps built with CUBISM for the 5.5 fim continuum, the 6.2 (im 
PAH feature, the 11.3/im PAH complex (11.2 + 11.3/jm PAH features), and the [Nc n] 12.81 fim emission line. For the PAHs and the 
emission line the nearby continuum was fitted and subtracted. The units of the maps are MJy sr — 1 for the 5.5 fim continuum image, and 
10 -4 erg cm -2 s -1 sr _1 for the PAH features and the [Nell] line maps. The orientation of the images is north up, east to the left. In all 
the maps, the (0,0) point corresponds to the approximate location of the nuclear region of IC 694. The images are shown in a linear scale. 
Only pixels with relative errors A/// < 0.2 are displayed. 



sociated background observations. We downloaded the 
BCD processed by the SpitzerlRS pipeline (version S15.3 
for SH and SI 7. 2 for LH). After removing rogue pixels 
with IRSCLEAN, we extracted a spectrum from each 
BCD using SPICE. In doing so, SPICE integrates the 
flux from the entire slit and performs an aperture cor- 
rection, which is dependent on the wavelength. We as- 
sumed point source flux calibration for the extraction. 
The final step was to median combine all the spectra for 
each of the two nuclei and to correct for fringing using 
IRSFPJNGE. 

3. ANALYSIS 

3.1. Spectral Maps 

For the SL data cubes we used cubism to construct 
spectral maps of the 6.2 fim aromatic (hereafter PAH) 
feature and the 11.3 /xm PAH complex (11.2 + 11.3 fim 
PAHs), the [Neil] 12.81 /im emission line, as well as the 



continuum at 5.5 /im. The integrated line flux maps 
are computed by integrating the line flux over the user- 
defined emission line regions with an average continuum 
removed. The 5.5 fim continuum map was constructed 
as the mean flux over a bandpass covering the spectral 
range of 5.3 — 5.7 fim. The map of the 6.2 fim PAH fea- 
ture was constructed by integrating the feature in the 
6.1 — 6.6 fim range and fitting the continuum between 5.3 
and 6.9 fim, whereas the 11.3 fim PAH feature was inte- 
grated between 11.1 and 11.8 fim and the corresponding 
continuum was fitted between 10.8 and 12.1 fim. Finally 
the [Ne n] emission line was integrated between 12.8 and 
13.1 fim, and the continuum was fitted between 12.2 and 
13.4 fim. All these wavelengths are observed values. 

Since cubism does not fit or deblend emission lines, the 
SL [Nell] 12.81 fim map includes some contribution from 
the 12.7 fim PAH feature complex. The SL spectral maps 
were then trimmed to the galaxy field of view (FOV) and 
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Fig. 2. — Maps of the line ratios (not corrected for extinction) observed with the SL module. The orientation of the images is north up, 
east to the left. The images are shown in a linear scale. The overlaid contours are the 6.2 (im PAH feature emission (left panel) and the 
11.3 yum PAH complex emission (right panel). 



rotated to the usual orientation of north up, east to the 
left. As discussed in §2.1 cubism produces both emission 
line and uncertainty maps. These spectral maps, not cor- 
rected for extinction, are shown in Fig. 1. In addition, 
we constructed maps of the observed 6.2 /xm to 11.3 /xm 
PAH ratio and the 11.3 /jm PAH to [Ne n] 12.81 (im ratio 
(see Fig. 2). Before constructing the line ratio maps the 
individual images were clipped to those regions with rela- 
tive errors of A/// < 0.2 (or S/N> 2), where / and Af 
are the integrated line flux and associated uncertainty, 
respectively, as derived by CUBISM. The typical S/N of 
the bright nuclei in the SL maps are between 400 and 
700 per pixel. 

For the bright emission lines present in the SH data 
cubes ([Siv]10.51/im, [Ne n] 12.81 /im, [Ne ill] 15.56 /xm, 
[Sm]18.71 /xm, and the H 2 S(l) line at 17.04 itm), we de- 
veloped our own idl routines to produce spectral maps. 
Briefly, we first automatically extracted ID spectra in 
2pixel x 2pixel boxes, by moving pixel by pixel along 
columns and rows until the entire FOV was covered. 
Each emission line was fitted with a Gaussian, and the 
continuum was fitted from adjacent spectral regions. In 
the case of the [Nell] 12.81 /im line we used two Gaus- 
sians, the second one to account for the nearby 12.7 /xm 
PAH feature. The lines were only fitted when the peak 
of the emission line was > 3c above the local continuum, 
where a is the standard deviation of the local continuum. 
If this requirement was not met, the pixel or pixels were 
masked out in the final spectral map. Since the 11.3 /xm 
PAH feature is not well reproduced with a Gaussian func- 
tion, we used CUBISM instead, as described for the SL 
data cubes, to construct the spectral maps. In this case, 
the 11.3 /xm PAH feature was integrated between 11.2 
and 11.7 /xm, and the continuum was fitted between 10.7 
and 12.3 /xm. 

After rotation of the SH spectral maps, we con- 
structed mosaics with the two pointings (Table 1) ob- 
tained for the system as shown in Fig. 3. We also con- 
structed maps of the [S in] 18.71 /im/[Neil]12.81 xtm and 



[Nein]15.56/xm/[Neii]12.81 /xm line ratios. We imposed 
a further restriction to include a given pixel in the map 
that the line errors of each emission line were less than 
25%. The maps of the two line ratios are shown in Fig. 4. 

3.2. Extraction of the ID Spectra 

For the two short wavelength IRS modules we first used 
CUBISM with the smallest (2x2 pixel) possible extrac- 
tion apertures on the data cubes before rotation. We ex- 
tracted ID SH and SL spectra for the brightest sources of 
Arp 299, namely A, Bl, C, and C (see Fig. 5). The opti- 
cal nucleus B2 is known to be faint in the mid-IR (Keto 
et al. 1997; Soifer et al. 2001), so we expect its con- 
tribution to the mid-IR spectrum of Bl to be minimal. 
However, there are Hn regions around Bl (see Fig. 5) 
that contribute to the SH and SL mid-IR spectra. 

Additionally we extracted SH spectra for two bright 
H II complexes. The first is in the southern spiral arm 
of IC 694, about 7" southeast of the nucleus, and it is 
clearly seen in the [Ne in] 15.56 /xm, [S in] 18.71 /xm, and 
[S iv] 10.51 /xm emission line maps in Fig. 3. The loca- 
tion of this region corresponds approximately to source 
F identified in Ha by Garci'a-Marm et al. (2006) as one of 
the youngest Hn regions in IC 694, based on the equiv- 
alent width (EW) of Ha. This region was also detected 
in the radio H92a emission line by Zhao et al. (1997). 
The second complex is located in NGC 3690, about 3" 
northwest from B2. In Fig. 5 the approximate locations 
of the SH extraction apertures are shown superimposed 
on the .H5T/NICMOS continuum-subtracted Paa image 
of Alonso-Herrero et al. (2000, AAH00 hereafter). The 
SH spectra for a 4.5" x 4.5" aperture of all these sources 
are shown in Fig. 6 (top and middle panels). The LH 
spectra of the three brightest regions A, B1+B2, and 
C+C were extracted with a ~ 13.4" x 13.4" aperture. 
We used CUBISM to extract matching SH spectra so that 
the full ~ 9.9 — 37 /xm spectral range was covered. The 
staring mode LH spectra of A and B1+B2 are shown in 
the lower panel of Fig. 6. 



The extreme star formation activity of Arp 299 
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Fig. 3. — Mosaic of the SH module observed (not corrected for extinction) continuum-subtracted emission lines maps, as well as the 
11.2 + 11.3 /jm PAH complex. The individual maps for the two pointings (Table 1) were built by fitting simultaneously the nearby continuum 
and Gaussians to the emission lines. Only pixels where the lines are detected at the 3<r level or higher are displayed. The units of these 
maps are 10~ 4 erg cm -2 s _1 sr — 1 . The spectral map of the 11.3 /^m PAH feature was built with CUBISM as explained in §3.1, and is shown 
in the same units as the other SH spectral maps. The orientation of all the images is north up, east to the left. The images are shown in a 
linear scale. 



Similarly, LL spectra with ~ 10.4" x 10.4" apertures 
were obtained for the brightest sources of Arp 299 (see 
Fig. 5), together with matching SL spectra so full ~ 
5 — 38 fim low-resolution spectra of A, B1+B2, and C+C 
were obtained, although there is some overlap between 
the B1+B2, and the C+C extraction apertures. Finally, 
we extracted a ID low-resolution (SL+LL) spectra cov- 
ering approximately the whole system using an LL aper- 
ture (before rotation) of 6 x 8 pixel (~ 30.5" x 40.6"). All 



the SL+LL spectra are shown in Fig. 7. 

3.3. Spectral features in the low-resolution spectra 

The low-resolution spectra (Fig. 7) for the three main 
components of the system (A, B1+B2, and C+C) as 
well as the integrated spectrum were fitted with PAHFIT 
(Smith et al. 2007a). pahfit uses a simple, physically- 
motivated model that incorporates starlight, thermal 
dust continuum in a small number of fixed tempera- 
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Fig. 4. — Mosaics of line ratios (not corrected for extinction) observed with the SH module. The images are shown in a linear scale. The 
overlaid contours are the [S III] 18.71 fim line emission (left panel) and the [Ncin]15.56 /im line emission (right panel). 



ture bins, resolved dust features and feature blends, and 
prominent emission lines, plus fully-mixed or screen dust 
extinction, including the silicate absorption bands at 9.7 
and 18 /im. pahfit is designed to work with the low 
spectral resolution data, and provides the best results 
when the full 5 — 40 /im spectral range is used. 

In Table 2 we give the pahfit results for the brightest 
emission line ratios (corrected for extinction), for the EW 
of the 6.2 /im PAH feature, and for the optical depth of 
the 9.7 /im silicate feature from the two dust configura- 
tions. The mid-IR extinction law used by pahfit can be 
found in Smith et al. (20071), The optical depth can be 
converted to optical extinction (Ay) using the conversion 
A v /T Si = 16.6 from Rieke & Lebofsky (1985). The two 
dust geometries provide fits of similar quality (see Ta- 
ble 2), with the screen model having lower silicate depths 
(between 2 and 3) than the fully-mixed model. This be- 
havior is similar to that found by Smith et al. (2007a) for 
the SINGS galaxies. Unlike the majority of the SINGS 
galaxies, the Arp 299 spectra are fitted with high val- 
ues of extinction, going from Ay ~ 12 to Ay ~ 33 mag 
(screen model geometry) for the different components of 
the system (see Table 2). 



with the SH module through a slightly larger extraction 
aperture (3pixel x 3pixel ~ 6.8" x 6.8"). The fluxes are 
about a factor of 1.6 larger indicating that some flux is 
missing from the 4.5" x 4.5" aperture, although the Paa 
emission (Fig. 5), radio continuum and line nuclear emis- 
sions (Aalto et al. 1997; Zhao et al. 1997) appear to be 
coming from a region less than 5" in diameter. 

The observed line ratios and [Ne n] line fluxes measured 
from the SH+LH spectra through 13.4" x 13.4" apertures 
are given in Table 4 for the three main components of 
the system. The fluxes and line ratios for IC 694-A, and 
the sum of NGC-3690-B1+B2 and NGC-3690-C+C' are 
in good agreement with those reported by Verma et al. 
(2003) from ISO large-aperture spectroscopy. We also 
note that the line ratios given in Table 4 are consistent 
within the uncertainties with those from the SL+LL data 
(Table 2) if they are corrected for extinction as derived 
from pahfit and the dust screen model. 

4. OVERALL MORPHOLOGICAL PROPERTIES 

In this section we describe the overall morphology of 
the different mid-IR spectral features, as well as the gen- 
eral physical and excitation conditions in Arp 299. 



3.4. Spectral features in the SH and LH spectra 

The fluxes of the brightest emission lines in the SH 
and LH spectra were measured using splot within the 
iraf noao environment. For all the emission lines the 
continuum was estimated from regions adjacent to the 
line. Most of the lines are well separated in the SH mod- 
ule, except for the [Ne n] 12.81 /im line. In this case, to 
get a good estimate of the line fluxes and the local con- 
tinuum, we fitted simultaneously three Gaussians to the 
H 2 S(2) line at 12.28 /im, the 12.7 /im PAH complex and 
the [Neil] line. 

The observed (not corrected for extinction) small- 
aperture SH line ratios in the six bright components de- 
scribed in §3.2 are listed in Table 3. We also give in this 
table the [Ne n] 12.81 /im line fluxes and errors, because it 
is the line with the largest errors as it is blended with the 
broad 12.7 /im PAH feature. Additional uncertainties in 
the fluxes and line ratios are due to centering errors and 
possible aperture corrections. As an example of the latter 
issue, in the notes to this table we also give the [Ne n] and 
[Ne in] fluxes for the nuclear region A of IC 694 measured 



4.1. Mid-IR continuum and fine- structure lines 

Figures 1 and 3 show that there is bright mid-IR con- 
tinuum and line emission from the two nuclei of the sys- 
tem (A and Bl), as well as from the C+C complex in 
the overlap region between the two galaxies. There is 
also emission in the spiral arms of IC 694 (especially in 
the southeast spiral arm), as well as in regions west of 
Bl and B2 in NGC 3690. These results are in agree- 
ment with previous work using optical, near-IR, mid-IR, 
millimeter, and radio continuum and emission lines (e.g., 
Gehrz et al. 1983; Aalto et al. 1997; Zhao et al. 1997; 
Satyapal et al. 1999; Sugai et al. 1999; AAH00; Soifer et 
al. 2001; Gallais et al. 2004; Garcia-Marm et al. 2006; 
Imanishi & Nakanishi 2006). 

In particular, the morphology of the fine-structure 
mid-IR [Ne in] 15.56 /im and [S in] 18.71 /im emission lines 
as well as the relative brightnesses from source to source 
resemble that of the ffST/NICMOS Paa emission line 
(Fig. 5, and AAH00) and other near-IR hydrogen recom- 
bination lines (Sugai et al. 1999; Satyapal et al. 1999). 
The hydrogen recombination lines are tracers of the 
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Fig. 5. — ffST/NICMOS images of the continuum-subtracted Paa emission at 1.875 (im (top panel) and of the continuum at 1.6/jm 
(bottom panel) of Arp 299 (see AAH00). The small squares (upper panel) and large squares (bottom panel) are the approximate locations, 
orientations, and sizes of the extraction apertures for the SH (~ 4.5" X 4.5") and the low-resolution (LL and matched SL apertures, 
~ 10.2" X 10.2") spectra, respectively. The orientation of the images is north up, east to the left. The images are shown in a square root 
scale. We also mark the positions of the two bright 1.6 /im continuum sources Bl and B2 in the lower panel. Note that there is very little, 
if any, Pact emission arising from B2 (upper panel). 



young ionizing stellar populations. The [S iv] 10.51 /^m 
line shows an overall morphology that is similar to Paa, 
but it is more affected by extinction, especially toward 
the nuclear region of IC 694 as we shall see in §7.1. 

4.2. The PAH Emission 

The PAH emission shows some morphological differ- 
ences when compared to the mid-IR fine structure lines, 
as can be seen from Figs. 1 and 3. In addition to the 
aforementioned bright sources in this system, there is dif- 
fuse bright PAH emission at 6.2 and 11.3 /itm (Fig. 1) in 
the interface region connecting the two galaxies, south- 
east of C+C. There is relatively bright PAH diffuse 
emission in the interface region where there is very lit- 



tle Paa emission (Fig. 5), possibly indicating that some 
PAH emission is not associated with the ionizing young 
stellar populations. 

Theoretical models predict that the relative strengths 
of the different PAH features depend on the properties 
of the dust grains, the charging conditions, as well as the 
starlight intensity (Draine & Li 2001). In particular, neu- 
tral PAHs are expected to show large 11.3 /im to 6.2 /xm 
ratios, whereas ionized PAHs have smaller ratios. The 
map showing the 11.3 /urn to 6.2 /im PAH ratio (Fig. 2, 
left panel), shows a range 10 from approximately 0.2 to 

10 The method used to construct the SL spectral maps for the 
PAH features is equivalent to a spline fitting, as opposed to full 
fitting of the spectra (see §3.3). In general Smith et al. (2007a) find 
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TABLE 2 

Results from pahfit for SL+LL spectra 



Source 


rsi 


[Nclll]l5.6 M m 
[NeIIJ12.8/am 


[ArllljS.O/im 
[ArII|7.0pin 


[BlVJlO.S^m 
[SIII]18.7Mm 


|fc!ni]33.5 M m 
[SIIIJ18.7/am 


[SlV]l().5 M m 
[ArIII|9.0/un 


|Sin]l8.7Mm 

[NeII|12.8/am 


EW 
6.2 fim PAH 


PAH6.2Mm 
PAH11.3nm 


x 2 


Screen Model 


A 


2.04 


0.06 


1.46 


0.47 


0.54 


0.37 


0.27 


1.36 


1.5 


1535 


B1+B2 


0.91 


0.32 


3.00 


0.03 


0.99 


0.04 


0.57 


0.22 


1.8 


1012 


C+C 


0.78 


0.34 


1.56 


0.32 


1.16 


0.46 


0.40 


0.50 


1.6 


670 


Integrated 


1.09 


0.36 


1.57 


0.15 


1.35 


0.38 


0.48 


0.65 


1.5 


1418 


Mixed Model 


A 


8.75 


0.05 


1.65 


0.35 


0.25 


0.33 


0.39 


1.85 


1.2 


1139 


B1+B2 


2.42 


0.32 


3.68 


0.06 


0.86 


0.12 


0.62 


0.22 


1.6 


823 


C+C 


1.98 


0.34 


1.68 


0.35 


1.07 


0.49 


0.43 


0.50 


1.5 


627 


Integrated 


2.98 


0.34 


1.84 


0.20 


1.23 


0.45 


0.54 


0.65 


1.4 


1170 



Note. — The extraction apertures are 10.4" x 10.4" for A, B1+B2, and C+C, and 30.5" x 40.6" for the integrated spectrum. 
The line ratios are corrected for extinction as derived from the measured 9.7 /tm optical depth and using the extinction law of 
Smith et al. (2007a). The EW of the 6.2 /im PAH feature, defined as EW= /(6. 2 /im PAH)//(6.2 /im continuum), are measured 
in /im. The value of x 2 given in the last column indicates the goodness of the fits to the spectra done with pahfit. 



0.9. From this figure it appears, as predicted by the mod- 
els, that the lowest PAH ratios are associated with the 
regions of high ionization (i.e., high [Ne ill] /[Neil] line 
ratios) in the galaxy. In comparison, the regions with 
high 11.3 /im to 6.2 /im PAH ratios are to the northwest 
of A and to the east of Bl where there is no evidence for 
bright Paa emission associated with ionizing stellar pop- 
ulations (Fig. 5). We note however, that both this ratio 
and the 11.3 /zm PAH/ [Neil] 12.81 ztm ratio are subject 
to extinction effects. This is because the 9.7 /zm absorp- 
tion feature is broad, and the 11.3 /zm feature is inside 
it (see figure 3 of Smith et al. 2007a for the mid-IR ex- 
tinction law of dust similar to that of the Milky Way) . 
Thus, caution should be exercised when interpreting ra- 
tios involving the 11.3 /im PAH feature (see discussions 
by Beirao et al. 2008 and Pereira-Santaella et al. 2009). 
For example, the observed values for the nuclear regions 
of both galaxies (A and Bl) are probably lower limits as 
both nuclei show significant values of extinction (§5 and 
§7.1). 

4.3. Molecular Gas 

The 17/zm H2 S(l) emission from IC 694 shows a 
more disk-like morphology compared to the mid-IR fine- 
structure emission lines, which mainly trace Hn emis- 
sion in the spiral arms. The disk scenario for IC 694 
is further confirmed by the observed H 2 velocity field 
(Pereira-Santaella et al. 2009). NGC 3690 shows H 2 
emission from the B1+B2 area, and the C+C complex. 
The H 2 S(l) emission of the whole system follows very 
closely the morphology of the molecular gas as traced by 
the 12 CO 1-0 line (Aalto et al. 1997; Casali et al. 1999), 
as also found by Sugai et al. (1999) for the near-IR H 2 
transitions. 

4.4. Physical and Excitation Conditions of the gas 

The [S ill] 18.71 /im/[Nell]12.81 /zm line ratio is almost 
insensitive to the hardness of the radiation field (i.e., age 

that for the 6.2 and 11.3 jum PAH features the full decomposition 
method measures intensities of about a factor of two higher than 
the spline method for both these two PAH features. We can then 
assume that the 11.3 /im to 6.2 (im PAH ratio map should give 
comparable results to the ratios fitted with pahfit. 



of the stellar population) and decreases with increasing 
densities, especially in the npi ~ 10 3 — 10 6 cm~ 3 range, 
depending on the ionization parameter (see Snijders et al. 
2007). The SH observed [S m] 18.71 /zm/[Ne 11] 12.81 /zm 
line ratio map (Fig. 4) shows spatial variations on phys- 
ical scales of ~ lkpc (FWHM ~ 2 pixels). Leaving aside 
the nuclear region of IC 694 11 , the observed ratio varies 
between approximately 0.2 and 0.7. 

Perhaps not surprisingly, those regions with the high- 
est [S ill] 18.71 /zm/ [Neil] 12.81 /im ratios do not coincide 
with the brightest mid-IR emitting regions. Rather they 
are located mostly in the southern part of IC 694, and 
to the north of the C+C complex and south of Bl. If 
this ratio traces the density (assuming a constant ion- 
ization parameter, cf §5 and 7) these would be the least 
dense regions in the system. The nuclei of the galax- 
ies (Bl and A) present the lowest ratios (see also Ta- 
bles 3, 4, and 5), and in particular A, the nuclear region 
of IC 694, would be the densest region, of the order of 
riH ~ 10 3 — 10 4 cm~ 3 (for solar metallicity, see Snijders 
et al. 2007). The significant change of the gas properties 
between the spiral arms and the nuclear region of IC 694 
observed in the [S ill] 18.71 /zm/[Neii]12.81 /zm spectral 
map is consistent with the variations of the molecular 
line ratios (Aalto et al. 1997). The molecular proper- 
ties of A require a population of unusually dense and 
warm clouds (see also Zhao et al. 1997). The C+C star- 
forming complex in the overlap region of the two galaxies 
presents [S ill] 18.71 /zm/[Ne 11] 12.81 /tm line ratios inter- 
mediate between those of the dense nuclear regions and 
the extranuclear regions of low surface brightness. 

The [Nelll]15.56/zm/[Nell]12.81/zm line ratio is sen- 
sitive to the hardness of the radiation field and has 
been quantitatively investigated by a number of works 
(e.g., Thornley et al. 2000; Martin-Hernandez et 
al. 2002; Verma et al. 2003; Rigby & Rieke 
2004; Snijders et al. 2007). This ratio has ad- 
vantages over other mid-IR hardness-dependent line 
ratios, in that it is less affected by differential ex- 
tinction than [S iv] 10.5 /tm/[S ill] 18.71 /zm, and can be 

11 the extinction corrections to the [S III] 18.7 fim line (see extinc- 
tion curve in Smith et al. 2007a and table 2 in Farrah et al. 2007) 
are non-negligible for this source (Table 2, and see also §7.1) 
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Fig. 6. — IRS high spectral resolution ID spectra for selected 
sources in Arp 299. Top and Middle: SH spectra of six bright 
sources in Arp 299 extracted from the spectral mapping data cubes 
with apertures of ~ 4.5" X 4.5" (see Fig. 5). Bottom: Staring 
mode LH spectra of the nuclear regions of IC 694 and NGC 3690 
extracted as point sources. The locations of the brightest emission 
lines are marked. 



measured from the SH spectra, unlike for instance 
[Ar lll]9.0 ^m/[Ar n]7.0 /im. 

Figure 4 shows the SH map of the observed 
[Ne in] 15.56 /im/[Neii] 12.81 /im line ratio, ranging from 
0.2 to approximately 0.9. The two nuclei show the small- 
est ratios, whereas the C+C complex and the southern 




10 



20 



Rest— frame X(/iim) 

Fig. 7. — Low spectral-resolution (SL+LL) spectra of the two 
nuclei of Arp 299 (A and B1+B2), the C+C complex, as well as 
the integrated spectrum for the interacting system. The extraction 
apertures (shown in Fig. 5) are 10.4" X 10.4" for A, B1+B2, and 
C+C, and 30.5" X 40.6" for the integrated spectrum. See §3.2 for 
details. 

spiral arm of IC 694, which are known to host some of 
the youngest regions in the system, show intermediate 
values. The highest ratios, however appear to be con- 
nected with regions of low surface brightness (e.g., the 
region about 10" northwest of Bl). It is not clear if the 
ratios in these regions reflect the youth of the stellar pop- 
ulation, as they do not seem to be associated with bright 
Hn regions detected in Paa (see Fig. 5), or rather result 
from other processes. For instance, there is bright diffuse 
UV emission to the west of sources Bl and B2 (Meurer et 
al. 1995). Moreover, a similar situation was observed in 
M82 where the [Ne in] 15.56 /xm/[Neii] 12.81 ^m line ratio 
increases at increasing distances form the galactic plane 
(Beirao et al. 2008). 

5. THE NUCLEAR REGION OF NGC 3690: THE OBSCURED 
AGN (Bl) AND SURROUNDING HII REGIONS 

We start out the detailed analysis of Arp 299 with 
the nuclear region of NGC 3690. This area contains 
two bright 1.6 /im continuum sources, namely Bl and B2 
(Fig. 5, lower panel). B2 is the optical nucleus, while 
Bl becomes the dominant IR source at A > 2/im in 
NGC 3690. B2 does not appear to be a bright source 
of mid-IR emission (Soifer et al. 2001). Direct evidence 
for an AGN in Bl comes from hard X-ray observations 
(Delia Ceca et al. 2002; Zezas, Ward, & Murray 2003; 
Ballo et al. 2004) and optical spectroscopy (Garcfa- 
Marm et al. 2006). Given the large X-ray column density 
(Ah = 2.5 x 10 24 cuT 2 , Delia Ceca et al. 2002; Zezas et 
al. 2003) toward this source the estimates of the X-ray 
luminosity range from Lx = 1-1 x 10 41 ergs _1 in the 
2-10 keV band (Ballo et al. 2004 from XMM observa- 
tions), to L x = 3 - 7 x lO^ergs" 1 in the 0.1 - lOkeV 
band (Zezas et al. 2003 from Chandra data). Regardless 
of the uncertain absorption corrections for this source, 
both estimates put it in the low-luminosity AGN cate- 
gory. Apart from the bright near-IR Bl and B2 sources, 
there is star-formation around Bl detected both in Paa 
(see Fig. 5) and in the UV (Gallais et al. 2004 and Meurer 
et al. 1995), indicating that these regions of star forma- 
tion must be relatively unobscured. 
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TABLE 3 

Observed Line Ratios from the SH spectral mapping data extracted with a 4.5" x 4.5" aperture. 



Source 


/([NeII]12.8/tm) 


[lNolll]15.6n.m 


[SlV)l().5 M m 


[SllrjlS. Vm 


H 2 M(2) 




[NoIIJ12.8nm 


|SIII]18.7fim 


[NeIIJ12.8/im 


H 2 S(1) 


|NcII]12.8nm 


IC 694-A 


11.0 ±2.2 


0.15 


0.13 


0.16 


0.65 


0.07 


IC 694-F 


2.0 ±0.4 


0.55 


0.22 


0.56 


0.42 


0.14 


NGC 3690-Bl(+B2) 


4.6 ±0.9 


0.44 


0.25 


0.33 




0.11 


NGC 3690-C 


6.0 ±0.6 


0.60 


0.28 


0.49 


0.45 


0.07 


NGC 3690-C 


3.0 ±0.3 


0.48 


0.30 


0.40 


0.71 


0.12 


NGC 3690-HII-NWB2 


2.8 ±0.6 


0.64 


0.21 


0.55 







Note. — The [Neil] fluxes are in units of 10 13 erg cm 2 s \ For comparison the 6.8" x 6.8" [Nell] and [Nem] fluxes for 
IC 694-A are 1.6 x 10~ 12 and 2.9 x 10~ 13 erg cm s , respectively. The fluxes and line ratios are not corrected for extinction. 
No aperture corrections have been applied. 



TABLE 4 

Observed Line Ratios from the SH+LH spectral mapping data extracted with a 13.4" x 13.4" aperture. 



Source 


/([NeII]12.8/tm) 


[NelUJ15.6/Mi> 


lSIVJ10.5(zm 


ism]18.7/jm 


H 2 S(2) 


H 2 S(i) 


lSIIIJ33.5Mm 


lSini34.8A.rn 


|NcII]12.8nm 


|SIIIJ18.7/jm 


[NcIIJ12.8Mm 


H 2 S(1) 


[NeIIJ12.8ftm 


[srn)i8.7^m 


[SIII]33.5^m 


A 


30.5 ± 6.1 


0.27 


0.18 


0.28 


0.18 


0.10 


0.8 


2.3 


B1+B2 


12.0 ±2.4 


0.52 


0.22 


0.40 






1.2 


1.4 


C+C 


18.8 ± 1.9 


0.58 


0.18 


0.51 






1.0 


1.2 



Note. — The [Nell] fluxes are in units of 10 13 erg cm 2 s 1 . The fluxes and line ratios are not corrected for extinction. No aperture 
corrections have been applied. 

5.1. AGN high excitation mid-IR emission lines 

The most direct way to identify an AGN in the mid- 
IR is through high excitation emission lines (e.g., Genzel 
et al. 1998; Sturm et al. 2002; Melendez et al. 2008). 
These high excitation lines are unlikely to be produced 
by star-formation, but they are not always detected in 
relatively bright AGN (e.g., Weedman et al. 2005). 
Among the high excitation lines, the [Ne v]14.32 /im 
and [O iv]25.89 /xm lines are found to be the bright- 
est in nearby AGN. We note however, that extended 
[O iv]25.89 /im line emission has been detected in some 
starburst galaxies, indicating that it can also be asso- 
ciated with star formation (Lutz et al. 1998). The 
[Ne v]14.32/xm is not detected in the SH spectra of Bl 
and surrounding regions (see Fig. 6, top panel). From the 
SH spectrum we placed a limit to the [Ne v] 14.32 /im line 
of < 3 x 10~ 14 erg cm -2 s _1 , consistent with the upper 
limit derived by Verma et al. (2003). The [O iv]25.89 /xm 
line on the other hand is detected in the staring mode 
LH spectrum (bottom panel of Fig. 6), and we measured 
a flux of 2 x 10~ 13 erg cm~ 2 s _1 (using an aperture cor- 
rection for point source). 

It is illustrative to compare Bl in NGC 3690 with 
the ULIRG NGC 6240, for which there is evidence for 
the presence of a double obscured AGN from X-ray ob- 
servations (Komossa et al. 2003), radio observations 
(Gallimore & Beswick 2004), and the detection of both 
[Nev]14.32/tm and [O iv]25.89 /im (Armus et al. 2006). 
The ratios between the 0.2 — 10 keV (corrected for absorp- 
tion) flux, and the [Nev] 14.32 /im and [O iv]25.89 /tm 
line fluxes for NGC 6240 are approximately 20 and 4, 
respectively. Using the 0.1 — lOkeV flux from the Chan- 
dra high spatial resolution data for the Bl source of 
Zezas et al. (2003) and the above flux ratios, we pre- 
dict line fluxes of /([NeV]) = 1-2 x 10~ 14 ergcur 2 s" 1 
and /([OIV]) = 6 - 9 x 10~ 14 erg cm" 2 s" 1 . The pre- 
dicted [Nev] flux is in agreement with the upper limit 




Rest— frame A(/im) 

Fig. 8. — SL spectra for a 3.7" X 3.7" aperture of the three main 
sources of Arp 299 compared to the average starburst spectrum of 
Brandl et al. (2006). All the spectra are normalized at rest- frame 
7.7 ^m. 

estimated above, whereas the predicted [O iv] flux is be- 
low the observed values. Similarly, using the correla- 
tion between the luminosities of the 2 — lOkeV emis- 
sion and the [O iv] 25.89 /im line found by Melendez et 
al. (2008) for a sample of Seyfert 1 galaxies, the esti- 
mated [O iv]25.89 /im luminosity in Bl is slightly below 
our own measurement, but within the scatter of the em- 
pirical correlation of Melendez et al. (2008). 

5.2. Hot dust emission in Bl and the composite nature 
of the mid-IR spectrum 

Bl appears as the brightest mid-IR continuum source 
of the entire system out to A ~ 20 /im (Figs. 6 and 7, 
and also Soifer et al. 2001; Charmandaris, Stacey, & 
Gull 2002). Gallais et al. (2004) remarked that this 
component shows a hot dust continuum at 5 /im, and 
they interpreted it as evidence for the presence of an 
AGN. This hot dust component is even detected down 
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Rest— frame A(/im) 

Fig. 9. — Rest-frame 5 — 8 (im spectra (thick black lines) of 
the three main sources of Arp 299 observed with the SL1 module 
and extracted with 3.7" X 3.7" apertures. The observed spectra 
are normalized at 6 ^m. For each component we show as dot-dash 
red lines the fitted AGN continuum, as dotted blue lines the fitted 
starburst component using the Brandl et al. (2006) template, and 
as thin solid green lines the sum of the fitted AGN + starburst 
components. We used a method similar to that of Nardini et al. 
(2008). 

to 2/im (Shier et al. 1996; AAH00). 

The hot dust continuum can be seen in the SL spectra 
of Bl in Fig. 8, where the spectra of the three bright- 
est sources are normalized at the 7.7 /im PAH feature. 
To quantify the hot dust emission we used a diagnostic 
method similar to that of Nardini et al. (2008) to sep- 
arate the AGN and starburst components using mid-IR 
spectra. This method is based on the close similarity of 
the 5 — 8 /im spectra of high metallicity starbursts (see 
Brandl et al. 2006). Thus, any excess in the 5 — 8/tm 
spectral region can be attributed to continuum emission 
from hot dust. This hot dust component can be taken 
as an indicator of an AGN (see Nardini et al. 2008 for 
further details). For our estimate we used two different 
templates, the starburst template of Brandl et al. (2006) 
and the star-forming ULIRG template of Nardini et al. 
(2008). By using two different templates, one can esti- 
mate the uncertainties associated with this method. We 
found that the AGN contribution at 6 /im is ~ 80 — 90% 
(see Fig. 9) for the 3.7" x 3.7" spectrum. This exam- 
ple illustrates the detection of a low-luminosity and/or 
heavily obscured AGN via the presence of a strong hot 
dust continuum at relatively short mid-IR wavelengths 
(3 — 6 /im, see also Risaliti et al. 2006) . 

Given that the IRS apertures include the contribu- 
tion from Bl (the AGN) and the surrounding Hn re- 
gions (see Fig. 5), it is reasonable to assume that this 
region be identified as a composite source (AGN + SB) 
based on several mid-IR diagnostics. To measure the 
apparent strength of the 9.7 /tm silicate feature, defined 
as Ssi = In /obs(9.7/im)// cont (9.7/im), we used a method 
similar to that outlined by Spoon et al. (2007) for sources 



dominated by PAH emission, fitting the continuum (see 
Fig. 8) as a power law with pivot wavelengths at 5.5 and 
14 /im. Indeed, the measured EW of the 6.2 /tm PAH 
feature (~ 0.2 — 0.3 /im) and the strength of the silicate 
feature (Ssi = —0.80 ±0.1) place source Bl in the Spoon 
et al. composite region, their IB region. Similarly, using 
the [Siii]34.8/im/[Sni]33.5/im vs. EW(6.2/iin PAH) di- 
agram from Dale et al. (2006), or the [Nev]/[Neii] vs. 
EW(6.2 /tm PAH) diagram from Farrah et al. (2007), Bl 
would be classified as a composite system. 

6. THE STAR-FORMING C+C' COMPLEX IN THE 
OVERLAP REGION 

Gehrz et al. (1983) identified the C+C' complex in the 
overlap region of the two galaxies, as the most luminous 
starburst in Arp 299, based on the observed Ha lumi- 
nosity. This region has been found to be not only very 
luminous, but also one of the youngest star-forming com- 
plexes in the system (e.g., Satyapal et al. 1999; Sugai et 
al. 1999; AAH00; Gallais et al. 2004; Garcia-Marin et al. 
2006), probably as the result of the interaction process. 
High angular resolution imaging at different wavelengths 
resolved the C+C complex into a number of star clusters 
and Hn regions (Meurer et al. 1995; Aalto et al. 1997; 
Casoli et al. 1999; AAH00; Soifer et al. 2001). CI, the 
brightest near-IR cluster (see AAH00 and Fig. 10) in the 
C region, and C are the most intense Paa sources in this 
region (see Fig. 10). Both appear to coincide with the 
bright radio and mid-IR sources in the overlap region. 

The detection of the Hei 1.70 /im line in C (Vanzi, 
Alonso-Herrero, Rieke 1998; Sugai et al. 1999) indi- 
cates that the temperature of the ionizing stars is high 
(T > 40, 000 K, see Vanzi et al. 1996; Forster-Schreiber 
et al. 2001). The observed Hell. 70 /im/BrlO line ratio 
(Sugai et al. 1999) would imply an age younger than ap- 
proximately 5 Myr for a solar metallicity instantaneous 
burst (see models by Rigby & Rieke 2004 and Snijders et. 
al 2007) or 7 Myr for super-solar models. The weak near- 
IR CO features (Ridgway, Wynn- Williams, & Becklin 
1994; Sugai et al. 1999; AAH00; Imanishi & Nakanishi 
2006) may indicate a lack of red supergiants confirming 
that this region is relatively young, but it is possible that 
the CO features are diluted by dust emission (see §6.3). 

6.1. Extinction and Physical Conditions 

The C+C complex appears to undergo varying degrees 
of obscuration based on the dust geometry as seen from 
the optical/near-IR color map of AAH00 (their figure 9). 
The dust distribution appears to be rather patchy and 
clumpy, especially near CI. This complicated dust con- 
figuration probably contributes to the range of extinc- 
tions (from a few magnitudes up to Ay — 37 mag) de- 
rived for this source using different optical and near-IR 
hydrogen recombination lines as well as different sized 
apertures (Satyapal et al. 1999; Sugai et al. 1999; 
AAH00; Garci'a-Marin et al. 2006). 

The SL+LL spectrum of C+C fitted with pahfit 
with a dust screen model provides an optical depth of 
r si = 0.8 (A v ~ 13 mag) over a 10.2" x 10.2" region. 
Since the dust distribution of this region appears to be 
patchy, we used the smallest possible SL1+SL2 apertures 
(3.7" x 3.7" ~ 750 pc x 750 pc) to extract spectra of four 
different regions in and around the C+C complex, to 
better isolate the different emitting sources. Apart from 
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Fig. 10.— Top: Contours (linear scale) of the HST/NICMOS 
continuum-subtracted Pact image of the C+C complex region. 
The squares show the approximate locations of the SL1+SL2 
3.8" X 3.8" extraction apertures for four different regions. We also 
mark the positions of bright near-IR continuum star clusters in this 
region as well as the optical nucleus B2 (see AAH00 for details). 
Bottom: SL1+SL2 spectra of the regions shown in the top panel. 
We also show the 10.2" X 10.2" SL+LL spectra that was fitted with 
PAHFIT. All the spectra are normalized at rest-frame 14 /im, one of 
the pivot wavelengths used to measure the apparent depth of the 
9.7 fim feature (see §6.1). 

sources C and C, the other two are a region at the in- 
terface of the galaxies southeast of C, characterized by 
faint Paa emission, and an H n region complex located 
approximately 3" northwest of B2. This H n region com- 
plex appears to be slightly older and less obscured than C 
and C, based on optical data (see AAH00 and Garcfa- 
Marm ct al. 2006). The approximate locations of the 
four regions are shown in Fig. 10. 

We measured the apparent strength of the 9.7 /^m sili- 
cate feature from the four SL1+SL2 spectra as described 
in §5.2. We found S 3i = -0.44 ± 0.05, -0.72 ± 0.10, 
-0.54 ± 0.05, and -0.25 ± 0.05, for C, C, the interface 
region, and the H n region complex NW of B2, respec- 
tively. These translate into values of visual extinction 
of Ay = 7 ± lmag and Ay = 12 ± 2 mag, to C and 
C, respectively. These values might only be lower limits 
as the broad PAH features at 7.7, 8.6 and 11.3 ^m fill 
the absorption feature resulting in an underestimate of 
the true optical depth. For instance, using the SL+LL 
10.2" x 10.2" spectrum of the C+C complex we obtained 
Ssi = —0.60 ± 0.1, whereas the optical depth fitted with 
pahfit for the same region is Tsi = 0.8 (Table 2). In the 



case of C, Sugai et al. (1999) using near-IR lines in the 
Brackett series derived Ay ~ 8mag for a 4.6"-diameter 
region, consistent with our estimate. 

The [Sni]33.48^m/[Siii]18.71^m line ratio can be 
used as an electron density estimator, as it has a small 
dependence on other properties such as the metallicity, 
ionization parameter, temperature of the stars, and the 
age of the ionizing stellar population (see e.g., Alexan- 
der et al. 1999; Snijders, Kewley, & van der Werf 2007). 
This line ratio is also most sensitive to densities in the 
range of ~ 10 2 — 10 4 cm~ 3 . Tables 2 and 4 give this ra- 
tio as obtained from the low (corrected for extinction) 
and high (not corrected for extinction) spectral resolu- 
tion data, respectively. The physical sizes covered by the 
extraction apertures are of the order of 2 — 2.7 kpc. 

The [Slii]33.5/Ltm/[Sill]18.7^m line ratios (Tables 2 
and 4) of the C+C complex are similar to those mea- 
sured in other star-forming galaxies on kiloparsec scales 
(e.g., [Sm]33.48Atm/[Sm]18.71/xm = 1.2 + 0.3 for the 
SINGS galaxies, Dale et al. 2007) and consistent with 
densities of the order of n c ~ 300 — 500 cm -3 (see also 
Verma et al. 2003). 

6.2. Mid-IR Emission Lines and the Star- Formation 
Properties 

The number of ionizing photons can be estimated from 
the [Ne n] and [Ne in] line luminosities using expressions 
such as those given by Roussel et al. (2006) and Ho 
& Keto (2007). These relations depend mainly on the 
neon abundance, and the radiation temperature. The 
dependence for the latter is through the fractional abun- 
dance of singly and doubly ionized neon. We took 
Trad = 4 x 10 4 K. We also assumed that most of the neon 
in Arp 299 is only singly and doubly ionized, consistent 
with young stars being the dominant source of ionizing 
radiation. In contrast to Ho & Keto (2007), who used 
a fraction of 0.6 for the proportion of ionizing photons 
that are absorbed by the gas as opposed to absorbed by 
dust or escaping, we take this fraction to be equal to 
one, that is, /i on = 1 in their equation 12. The last piece 
of information needed to estimate the ionizing photons 
is the neon abundance. Using ISO mid-IR spectroscopy 
Verma et al. (2003) measured super-solar neon abun- 
dances in the two components of Arp 299. Further ev- 
idence for super-solar abundances in this system comes 
from the oxygen abundance of log(0/H) + 12 = 8.92 
derived from the optical spectroscopy of Moustakas & 
Kennicutt (2006) and using the Kobulnicky & Kewley 
(2004) calibration (John Moustakas, private communi- 
cation). This oxygen abundance is consistent with the 
range derived by Rupke, Veilleux, & Baker (2008) for a 
large sample of LIRGs in the local universe. We com- 
puted the number of ionizing photons using the expres- 
sion given by Ho & Keto (2007) but using the neon abun- 
dances derived by Verma et al. (2003) for NGC 3690, 
that is, [Ne/H]= 1.8 x 10~ 4 instead of the solar value of 
[Ne/H]= 9.9 x 10~ 5 . 

Under the assumptions stated above we find -/Vl vc ~ 
3 x 10 53 s- 1 and N hyc ~ 2x 10 53 s" 1 for C and C, respec- 
tively. The value for C is a factor of 10 lower than the es- 
timate from Paa using Ay = 15mag derived by AAH00. 
Using the values of Ay — 7mag and Ay = 12mag derived 
for C and C, and the Br7 fluxes of Sugai et al. (1999), 
which were measured through 4.6"-diameter apertures, 
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we get N hyc - 2 x lO^s" 1 and N Lyc ~ 1 x lO^s" 1 , 
respectively. The discrepancy between A^Lyc estimated 
from extinction-corrected Br7 and from the neon fine 
structure lines could be explained if some of the ionized 
gas is at densities exceeding n cr ;t ~4x 10 5 cm~ 3 for the 
neon lines. This situation would arise if ultra-compact 
H II regions play an important role in the star forma- 
tion in regions C and C, as suggested by Rigby & Rickc 
(2004) for luminous young starbursts in general. We re- 
turn to this possibility at the end of §6.3. 

The ionization parameter (defined as q = 
-^Lyc/(47ri? 2 rti on )) can then be estimated measuring ap- 
proximate sizes of the nebulae from the HST/KLCMOS 
Paa images: R <~ 300 pc and R <~ 250 pc, for C and 
C, respectively, and taking rin = 300 cm~ 3 . We find 
g ~ 0.8 - 4 x 10 8 cms _1 for both C and C, using the 
neon line and hydrogen recombination line estimates of 
iV Lyc - 

We now make the simplifying assumption that the ion- 
ized gas can be divided entirely between high and low 
density regions. We use the mid-infrared fine structure 
lines to diagnose conditions in the low-density region. 
Diagrams involving mid-IR line ratios of the same ele- 
ment can be used to estimate the age of stellar popu- 
lation, although there is some degeneracy between the 
age and the ionization parameter. We use the Snijders 
et al. (2007) predictions of the age evolution of mid- 
IR fine-structure emission lines. These are calculated by 
combining the outputs of Starburst99 (Leitherer et al. 
1999) and the photoionization code Mappings (Dopita 
et al. 2000). The stellar populations are assumed to 
have been formed in an instantaneous burst of star for- 
mation with a Salpeter IMF (between 0.1 and 100 M Q ). 
Given the abundances measured in Arp 299 we used the 
2Zq models. 

Figure 11 is a [Ne in] 15.56 ^m/[Neii]12.81 ^m vs. 
[S iv] 10.51 /zm/[S in] 18.71 /im line ratio diagram with the 
model predictions of Snijders et al. (2007). We show 
the models for two nebular densities, rin = 10 2 cm~ 3 , 
the low-density limit typical of star-forming galaxies (see 
above), and «h = 10 4 cm -3 to account for the high den- 
sity nuclear region of IC 694 (see §7.1). We also show 
the model predictions for three ionization parameters 
representative of those measured in different sources in 
Arp 299 (see also §7.2 for A). 

The line ratios (measured with the 4.5" x 4.5" aperture) 
for C and C together with those of other sources 12 in the 
system are plotted in Fig. 11. The right panel gives the 
time evolution for the first few million years. The line 
ratios drop very fast as the most massive stars move off 
the main sequence with a minimum at about 3 Myr. Af- 
ter that the first Wolf-Rayet (W-R) stars begin to appear 
the ratios increase again (left panel of Fig 11) for the du- 
ration of the W-R phase (out to approximately 7 Myr for 
the super-solar models, Leitherer et al. 1999; Snijders et 
al. 2007). Given the measured ionization parameters for 
C and C', we can rule out the youngest ages for these 
two sources. The mid-IR line ratios in Fig. lis could be 

12 We also plotted in Fig. 11 the line ratios of Bl. We caution 
the reader that with the spatial resolution of the IRS data we 
cannot estimate what fraction of the emission line fluxes originates 
from the AGN and what fraction from the surrounding H II regions. 
Thus we do not discuss the star- formation properties of Bl based 
on this diagram 



reproduced with ages 4 — 7 Myr, using the 2Zq models. 
A summary of the mid-IR derived properties of the low 
density zones in C and C' sources is given in Table 5. 
It is interesting to note however, that given the derived 
ranges of q and ages, the [S in] 18.71 /xm/[Neii] 12.81 /im 
line ratios can only be reproduced with relatively high 
nebular abundances (njr > 10 4 cm~ 3 ). 

6.3. W-R features, hot dust and radio emission 

To further explore the possible presence of young stel- 
lar populations in this region, we looked for the W-R 
broad He II feature at 4686A in the optical integral field 
spectroscopy (IFS) of Garcia-Marm et al. (2006). The 
detection of this feature provides direct evidence for the 
presence of a large population of hot, massive and young 
stars (Vacca & Conti 1992). The optical spectra of a 
few bright sources in Arp 299 are shown in Fig. 12. The 
W-R broad He n feature at 4686A is clearly detected in 
C and C', in the Hn regions surrounding Bl, and tenta- 
tively detected in the H II region complex NW of B2. The 
high metallicity models predict more W-R stars over an 
extended lifetime compared with solar metallicity (Lei- 
therer et al. 1999). Specifically the 2Zq models predict 
that stars enter the W-R phase at about 2.5 Myr and 
dominate the radiation field out to 7 Myr (see Snijders 
et al. 2007 for more details). 

C and C' are also bright radio sources (Gehrz et al. 
1983; Neff et al. 2004) with observed radio spectral in- 
dexes 13 of a = —0.65 and a = —0.50, respectively (Neff 
et al. 2004). Models for the radio thermal emission of 
young stellar populations (< 6 Myr) and solar metallicity 
predict a > —0.1 (Perez-Olea & Colina 1995), indicating 
that some fraction of the radio emission must originate 
from supernovae. Using the N-^ yc derived in §6.2 from 
Br7, the predicted thermal emission (see Colina & Perez- 
Olea 1992) would account for approximately 30% of the 
observed emission at 1.4 GHz (Gehrz et al. 1983) from 
C through a 4" x 4" aperture. We cannot make mean- 
ingful comparisons at the other radio frequencies (4.9 
and 8.4 GHz), because the measurements of Gehrz et al. 
(1983) were with beam sizes (~ 1") much smaller than 
our typical angular resolution (see §2.1). Given this, it 
is possible that both C and C are just now beginning a 
significant supernova activity, consistent with the lack of 
young radio supernovae (Neff et al. 2004). All these facts 
confirm an age of C and C' of approximately 4— 7 Myr de- 
rived from the mid-IR line ratios, rather than extremely 
young bursts (see §6.2). 

Adding to this rather complicated situation, is the 
presence of a hot dust component, which accounts for 
approximately 30 — 40% of the 6 /xm emission in C within 
the 3.7" x 3.7" SL aperture (see Fig. 9). A hot dust con- 
tinuum has also been found based on the observed EW 
of the 3.3^m PAH feature (Imanishi & Nakanishi 2006). 
Hot dust is uncommon in starbursts of the luminosity of 
C (Ridgway et al. 1994). However, it is unlikely that 
this dust is associated with a buried AGN, as there is no 
other evidence for the presence of an active nucleus. 

The hot dust can be associated with ultra-compact re- 
gions, whose presence is also indicated by the discrepancy 
in A^Lyc deduced from Paa and the neon fine structure 

13 defined as S„ oc v a , where S v arc the flux densities at 4.9 and 
8.4 GHz 
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Fig. 11. — Mid-IR line ratio diagrams. The symbols are as follows, the filled triangle is the extinction-corrected IC 694-A (using the 
scaling factors in table 2 of Farrah et al. 2007 for the Li & Draine 2001 extinction law); the star symbol is NGC 3690-B1; filled circles are 
the C and C sources (marked for clarity in the upper right panel); the open circles are the F region in IC 694, and the Hn complex NW 
of B2 in NGC 3690. All the data are from the 4.5" X 4.5" SH spectra. The solid lines are the Snijders et al. (2007) models for 2Zq and 
different ionization parameters. The top panels are the njj = 10 2 cm -3 models and lower panels are the nu = 10 4 cm -3 models. The right 
panels show the model evolution for 0, 1, and 2Myr, whereas the left panels are the model evolution from 3.5, 4, 5, 6, and 7Myr. The 
dotted lines connecting the different models represent the "isochrones" in these diagrams. 



lines. It is likely that what we are calling source C is 
in reality a collection of different sources with different 
ages, as already hinted by the detection of a number of 
near-IR continuum sources in the large Paa nebula (see 
Figs. 5 and 10). Source C is reminiscent of the nuclear 
region of the W-R galaxy Henize 2-10, which hosts young 
clusters, ultra-dense Hn regions, and possibly supernova 
remnants within the central few hundred parsecs (see Ca- 
banac, Vanzi, & Sauvage 2005; Martin-Hernandez et al. 
2006). 

The discovery that ultra-compact regions play an im- 
portant role in the young star-forming regions in C and 
C supports the hypothesis of Rigby & Rieke (2004). 
They proposed that large numbers of such regions would 
account for the lack of high-excitation emission lines in 
starbursts without imposing restrictions on the initial 
mass functions in them. 

7. THE NUCLEAR REGION OF IC 694 (A) 



From the early studies of this interacting system, the 
star-formation properties of the nuclear region of IC 694 
were found difficult to model, in part due to the strong 
dust extinction that obscures most of the Ha emission 
and in part due to the presence of bright radio emission 
(Gehrz et al. 1983). It also became apparent that this 
source alone might be responsible for about half of the 
total IR luminosity of the system (Gehrz et al. 1983; 
Joy et al. 1989; Keto et al. 1997; Charmandaris et al. 
2002). Near-IR imaging and spectroscopic observations 
were able to partially overcome the strong obscuration 
to the nucleus of IC 694. Sugai et al. (1999) and AAH00 
showed that this nuclear region is affected by strong ob- 
scuration, and that the observed properties are the re- 
sult of an extended period of star formation (see also 
Shier et al. 1996). In fact, the star-formation activity 
of source A appears to be older than that of the other 
bright components of the system. This conclusion was 
later reinforced by the detection of a number of radio 
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Fig. 12. — Optical spectrum of a number of sources in Arp 299 
from the IFS of Garcfa-Marin et al. (2006) extracted with their 
large fiber of 3" (~ 610 pc) in diameter to match as much as pos- 
sible the small Spitzer apertures. The W-R broad He II feature at 
4686A is clearly detected in C and C, and marginally detected in 
the H II region complex NW of B2. 

sources in the central 0.6" of A, which could only be 
explained as supernova remnants (Neff et al. 2004; Ul- 
vcstad, Johnson, & Neff 2007) from a relatively evolved 
star formation (6 — 8 Myr after the peak of the activity) . 
Other works suggested the presence of an obscured AGN 
in the nucleus of IC 694. This possibility is discussed in 
§7.3. 

7.1. Extinction and Physical Conditions 

The overall shape of the low-resolution spectra of the 
nuclear region of IC 694 appears to be remarkably dif- 
ferent from the other components of Arp 299 (see Fig. 7 
and 8, and Gehrz et al. (1983). The mid-IR spectrum 
of A shows a much deeper absorption at 9.7 /im (see also 
Gallais et al. 2004), and 18 /im than the other compo- 
nents in the system and than the average spectrum of 
starburst galaxies (Brandl et al. 2006). The depth of 
the silicate feature is comparable to the values observed 
in local ULIRGs (see e.g., Spoon et al. 2007; Armus 
et al. 2007; Farrah et al. 2007). pahfit fits a 9.7 /xm 
apparent optical depth of r s j = 2.04 for the dust screen 
model. Using the Rieke & Lebofsky (1986) extinction law 
this optical depth implies an extinction of Ay ~ 34 mag 
(see also Gallais et al. 2004), comparable to, although 
higher than, the near-IR estimates (e.g., Beck, Turner, 
& Ho 1986; Sugai et al. 1999; AAH00). We note, how- 
ever, that the observed apparent optical depth of A and 
other ULIRGs cannot be reproduced with a foreground 
dust screen, and they require the dust to be in a deeply 
embedded and optically thick smooth distribution (Lev- 
enson et al. 2007). 

The nuclear region of IC 694 over a 2 kpc scale shows 
a much lower [Slii]33.48/zm/[Slll]18.71 //m line ratio 
(corrected for extinction) than the nuclear region of 
NGC 3690 (Bl) and the C+C complex. Using the pho- 
toionization models of Snijders et al. (2007), the ob- 
served sulfur ratios of the central 2 — 2.7 kpc of IC 694 
would imply densities of the order of n c ~ 1— 5xl0 3 cm~ 3 
(see also models of Alexander et al. 1999). This indicates 
that the nuclear region of IC 694 has a density higher 
than the typical values inferred for starburst galaxies 



from mid-IR spectra (n e ~ 300 — 400 cm 3 , Verma et 
al. 2003; Dale et al. 2006). 

The high density of the nuclear region of IC 694 derived 
from the [S in] lines is in agreement with the modelling 
of the radio H92a transition of Zhao et al. (1997). They 
found values of n c ~ 5 x 10 2 - 5 x 10 4 cuT 3 at the very 
center and n c ~ 5 x 10 3 — 10 4 cm -3 in the nuclear exten- 
sion about 2" southeast of the center (also seen in the 
HST/mCMOS Paa map, Fig. 5). Aalto et al. (1997) 
also found that the high 12 CO/ 13 CO molecular line ratio 
of the A nuclear region implied a population of unusually 
dense and warm clouds, due primary to small or moder- 
ate optical depth of the 1-0 12 CO transition. 

7.2. Mid-IR Emission Lines and the Star- Formation 
Properties 

The number of ionizing photons arising from the cen- 
tral region of IC 694 has been estimated by a number 
of IR and radio methods, and is in the range N^ yc ~ 
1 - 4 x lO^s" 1 (Shier, Rieke, & Rieke 1996; Zhao et 
al. 1997; Sugai et al. 1999; AAH00). We can obtain 
an independent estimate from the extinction-corrected 14 
[Ne n] 12.81 /jm and [Ne in] 15.56 /jm luminosities mea- 
sured with the SH 6.8" x 6.8" aperture, as explained in 
§6.2 and using the neon abundance inferred by Verma et 
al. (2003). We find A Lyc - 0.8 x lO^s" 1 from the neon 
lines. As found for C and C the number of ionizing pho- 
tons from the neon lines is a factor of a few lower than 
other estimates from hydrogen lines. The corresponding 
range of ionization parameters is q ~ 1 — 5 x 10 s cms -1 
using a size of the nebula of R ~ 200 pc measured from 
the #ST/NICMOS Paa image), n H = 1000cm" 3 , and 
the number of ionizing photons from the neon lines and 
the Br7 measurements of Sugai et al. (1999). 

The observed mid-IR line ratios 
([Neni]15.56 /zm/[Neii]12.81 ^m and 
[Slv]10.51/im/[Sin]18.71 /im) of source A are shown 
in Fig. 11 with the filled triangle symbol. Since the 
highest possible ionization parameter for A would be 
q ~ 4 x 10 cms (i.e., using the number of ionizing 
photons from the Br7 line of Sugai et al. 1999), we 
can rule out the youngest age for A (1 — 2 Myr, lower 
right panel of Fig. 11). The mid-IR line ratios from A 
would be better explained with an evolved starburst 
> 6 — 7 Myr (lower left panel of Fig. 11). This is con- 
sistent with the age derived from the radio and near-IR 
properties. We also note that models with constant or 
quasi-constant (e.g., exponentially decaying with long 
timescales) star formation can also be excluded for this 
source as such models would predict high mid-IR line 
ratios for all ages (see Thornley et al. 2000; Rigby &; 
Rieke 2004). 

7.3. No mid-IR evidence for an obscured AGN 

Early suggestions for the presence of an obscured AGN 
in A were based on the compact nature of the radio emis- 
sion (e.g., Sargent & Scoville 1991; Lonsdale, Lonsdale, 
& Smith 1992). Ballo et al. (2004) measured an X-ray 
luminosity for A similar to that of Bl (see §5), and found 
it could not be accounted for by star formation alone. In 
contrast, Zezas et al. (2003) were able to explain this 

14 Ay = 34 mag and the Rieke & Lebofsky (1985) extinction law 
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TABLE 5 

Summary of Mid-IR derived properties of the main components of Arp 299 



Source 


A v 


n c 


Age 


A^Lyc 




<? 


Hot dust 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


A 


34 


1000 - 5000 


> 7 


8 


1 - 


- 5 x 10 s 


10-20% 


Bl 


14 


250 - 400 










80-90% 


C 


7 


250 - 400 


4-7 


3 


0.8 


- 4 x 10 8 


30-40% 


C 


12 


250 - 400 


4-7 


2 


0.8 


- 4 x 10 8 





Notes. Column (1): Source. Column (2): Optical extinction in magnitudes as derived from the optical depth of the 9.7 (im absorption 

feature. For A the value is that obtained from pahfit, whereas for the other sources the apparent depth was measured as explained in 

§5.2. Column (3): Electron density in cm -3 from the [S m]33.5 /zm/[S nr]18.7/im line ratio. Column (4): Age estimate in Myr from 
mid-IR line ratios. Column (5): Number of ionizing photons in units of 10 53 s — 1 estimated from the [Nell] 12.81 fim and [Ne III] 15.56 fim 
emission lines. The two neon line fluxes were measured in 4.5" X 4.5" apertures for C and C', and for 6.8" X 6.8" for source A. The neon 

abundances used are 2.2 and 1.7 solar for A, and C+C, respectively (Verma et al. 2003). Column (6): Range of possible ionization 
parameters in units of cm s _1 , using the two estimates of the number of ionizing photons, from neon lines and Br7 measurements of Sugai 
et al. (1999). Column (7): Hot dust contribution to the 6 /an luminosity using a method similar to that described by Nardini et al. (2008). 
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Fig. 13. — Comparison of the low spectral-resolution spectra 
of the integrated emission of Arp 299 with the average starburst 
template of Brandl et al. (2006), the average spectrum of PAH- 
dominated local ULIRGs (the 2C class of Spoon et al. 2007), and 
the average spectrum and corresponding la dispersion of high- 
z ULIRGs from the Farrah et al. (2008) sample. The IR (1 - 
1000 /im) luminosity range of the galaxies included in this high-z 
template is ~ 8 X 10 12 to 6 X 10 13 Lq. The spectra have been 
scaled to match the peak of the 6.2 fim PAH feature. 



luminosity in terms of high mass X-ray binaries as esti- 
mated from the number of ionizing photons. The water 
megamaser properties of this nucleus do not exclude the 
possibility of an AGN in A (Tarchi et al. 2006). 

As for Bl we did not detect the high excitation 
[Nev] 14.32 fim emission line, typical of AGN. The 
[Oiv]25.89 /urn is detected in the LH staring mode spec- 
trum (bottom panel of Fig. 6) with a measured flux of 
1 x 10 _14 ergcm -2 s _1 . Unfortunately the limited spatial 
resolution of the LH data does not allow us to determine 
whether the integrated [O iv] emission overr the LH slit 
(~ 11" x 22") is extended or not. Using the rest-frame 
5 — 8 /im spectrum and the method described in §5.2, 
the contribution from hot dust at 6 /im is estimated to 
be approximately 20%. Given the uncertainties of the 
method, this contribution is formally compatible with 
pure star formation (see Fig. 9). Finally the location of 
A on the diagram of EW of the 6.2 /im PAH feature vs 
the depth of the silicate feature (Spoon et al. 2007) is 
similar to that of M82, and thus consistent with obscured 
star formation. 



8. ARP 299 IN THE CONTEXT OF LOCAL AND HIGH-Z 
IR-BRIGHT GALAXIES 

The unprecedented sensitivity of IRS on Spitzer is mak- 
ing it possible to characterize the mid-IR spectra of IR 
bright-galaxies at z > 1. Perhaps one of the most sur- 
prising results is that a large fraction of these high-z 
IR-selected galaxies are classified as ULIRGs in terms of 
their IR luminosities, but their mid-IR spectra are more 
similar to those of local star- forming galaxies (e.g., Yan et 
al. 2007; Farrah et al. 2008; Rigby et al. 2008). That is, 
ULIRGs at high-z show strong PAH features and mod- 
erate depths of the 9.7 /im silicate feature. This is in 
contrast with their local counterparts that tend to show 
very deep silicate features and moderate to low equiva- 
lent widths of the PAHs (e.g., Spoon et al. 2006; Armus 
et al. 2007; Farrah et al. 2007). 

Figure 13 shows a comparison of the integrated spec- 
trum of Arp 299 with the average template of starbursts 
from Brandl et al. (2006). It is immediately clear from 
this figure that the two mid-IR 5 — 14 /im spectra are 
quite similar, at least in terms of the PAH features, de- 
spite the different range of IR luminosities of the galaxies 
included in the starburst template 15 . It is also clear that 
the integrated mid-IR spectrum of Arp299 is not dom- 
inated by the strongly obscured source in the nuclear 
region of IC 694 (see Fig. 7). The apparent strength of 
the 9.7 /im silicate feature of the integrated spectrum is 
Ssi = —0.79 ± 0.05, similar to the typical values shown 
by starburst galaxies (Spoon et al. 2007; Brandl et al. 
2006), whereas for source A we find S Si = -2.03 ± 0.10. 
We also show in this comparison the average template of 
local PAH-dominated ULIRGs (the 2C class of Spoon et 
al. 2007), which is similar in terms of the PAH features 
but it has a deeper silicate feature. The other domi- 
nant class of local ULIRGs (not shown in Fig. 13), the 
3C class in the Spoon et al. (2007) classification, has a 
much deeper silicate feature and low EW of the PAH fea- 
tures, and it is even more different from the integrated 
spectrum of Arp 299 . 

We also show in Fig. 13 the average spectrum of the 
sample of high-z ULIRGs (and even higher IR luminosi- 
ties) of Farrah et al. (2008). The galaxies in this sam- 
ple were selected such that they showed a strong 1.6 /xm 
feature and bright 24 /im fluxes, and are mostly in a 

15 This template is the average of 13 starburst galaxies, of which 
two galaxies are just above the LIRG limit, whereas the rest of the 
galaxies have IR luminosities between 6 X 10 9 Lq and 6 X 10 10 Lq. 
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narrow rcdshift range around z = 1.71. The similarity 
between the integrated spectrum of Arp 299 and these 
high-z ULIRGs is quite remarkable. One of the favored 
explanations for the observed mid-IR spectra of high-z 
ULIRGs is that star formation is diffuse and extended, or 
distributed in multiple dusty star-forming regions spread 
over several kiloparsecs. Arp 299 behaves in this manner 
and may be a local counterpart of the star formation pro- 
cesses taking place in high-z IR-bright galaxies (see also 
Charmandaris, Le Floch, & Mirabel 2004), albeit with a 
lower IR luminosity and possibly higher metallicity. This 
is in contrast with local ULIRGs where most of the star- 
formation is occurring in the nuclei of the galaxies on 
sub-kiloparsec scales (Soifer et al. 2000). 

9. CONCLUSIONS 

We have presented Spitzer/IKS spectral mapping of the 
central - 45" - 9 kpc of Arp 299 (IC 694 + NGC 3690). 
The IRS data combined with #5T/NICMOS imaging 
and optical ground-based IFS allowed us to study in 
detail the morphological distribution of mid-IR spectral 
features as well as the mid-IR properties of the bright- 
est sources in the system. We used a number of mid-IR 
indicators to derive the star formation properties (ages, 
number of ionization photons, and ionization parame- 
ters) and physical conditions of selected regions in the 
system. The main conclusions of this study are: 

• The spatial distribution of the bright mid-IR emis- 
sion lines follows well the ionizing stellar popula- 
tions as probed by the ffST/NICMOS Paa emis- 
sion. High [Ne in] 15.56 /im/[Nen] 12.81 /um and 
[Siv] 10.51 ^m/[S in] 18.71 /jm line ratios mark the 
location of the youngest regions, that is, the C+C' 
complex and H n regions in the disks of the galax- 
ies. 

• Most regions in the system show 
[Sin]33.48/xm/[Slll]18.71/im line ratios similar to 
those typical of star-forming galaxies, implying 
densities of the order of n c ~ 250 — 500 cm~ 3 . 
Only the ratios of the nuclear region of IC 694 
indicate higher densities, n c ~ 1000 — 5000 cm~ 3 , 
in agreement with previous results for this source 
from radio hydrogen emission and molecular CO 
line ratios. 

• The extinction as estimated from the optical depth 
of the 9.7 fim silicate feature is found to vary signifi- 
cantly from region to region in Arp 299. For the nu- 
clear region of IC 694 we measured rsi = 2 or Ay = 
34 mag, comparable to the optical depths measured 
in local ULIRGs. The other bright sources of the 
system (Bl, C, and C) show more moderate values 
of the extinction ranging from 7 to 14 mag. 

• Both 6.2 /xm and 11.3 /jm PAH emission is detected 
in the bright sources of the system as well as at the 
interface region connecting the two galaxies. There 
is a tendency for the observed 11.3 to 6.2 ^m PAH 
feature ratio to be higher in regions of low ion- 
ization, although the observed 11.3 ^m PAH fea- 
ture emission can be significantly reduced in re- 
gions with elevated extinction. 



• The mid-IR line ratios suggest that youngest stel- 
lar populations (~ 4 — 7Myr, using 2Zq models) 
are mainly located in the overlap region of the two 
galaxies (the C+C complex), the southern spiral 
arm of IC 694, and the H n region complex to the 
northwest of the nuclear region of NGC 3690. The 
detection of optical W-R features in C and C pro- 
vides further evidence for the youth of these re- 
gions. The extended lifetimes of super-solar metal- 
licity W-R stars means that they can dominate the 
radiation field to older ages (out to ~ 7Myr). In 
contrast, the mid-IR properties of the nuclear re- 
gion of IC 694, in particular the low [Nelll]/[Neii] 
line ratios, confirm that the star formation activity 
there is slightly older and more extended in time 
than in the aforementioned regions. 

• Evidence for hot dust emission in the 5 — 8 /jm spec- 
tral region for Bl, the nuclear region of NGC 3690. 
Bl hosts an X-ray and optically identified low- 
luminosity AGN. The AGN is likely to be respon- 
sible for ~ 80 — 90% of the 6 /xm emission (within 
a 3.7" x 3.7" aperture). 

• In C, the overlap region between the two galaxies 
and a site of very recent star formation, we found 
both hot dust and also evidence that a portion 
of the ionized gas may be in very dense regions, 
above the critical densities for [Neil] and [Nem]. 
The hot dust continuum contributes <~ 30 — 40% at 
6 /an. Given that there is no other indication for 
the presence of an AGN in C, we interpreted it as 
produced by dust heated by a young star cluster 
or clusters. Both of these results may be consis- 
tent with the presence of a substantial population 
of ultra-compact Hn regions. 

• The integrated mid-IR spectrum of Arp 299 is sim- 
ilar to that of local starbursts despite its strongly 
interacting nature and high IR luminosity. This is 
because the star formation in this system is spread 
across at least 6-8 kpc, with a large fraction taking 
place in regions of moderate mid-IR optical depths. 
It is only the nuclear region of IC 694 that shows 
the typical mid-IR characteristics of ULIRGs, that 
is, very compact and embedded star formation re- 
sulting in a deep silicate feature and moderate EW 
of the PAHs. We also find that the integrated 
5 — 14 /im spectrum of Arp 299 is similar to the 
average spectrum of the high-z sample of Farrah et 
al. (2008). This suggests that the Arp 299 system 
may represent a local example, albeit with lower IR 
luminosity and possibly higher metallicity, of the 
processes occurring in high-z IR-bright galaxies. 
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